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During  the  final  period  November  1981  -  September  1983, 
the  contract  continued  to  support  Dr.  U.  Spring  as  a  senior 
research  associate,  followed  by  a  gap  of  3  months,  then 
supported  Dr.  H.  T.  Williams  as  a  senior  research  associate. 

This  report  is  presented  in  two  parts,  I  and  II.  Part  I 
has  now  been  published  in  Cold  Regions  Science  and  Technology, 
6^  (  1983)  185-193,  being  a  continuation  of  my  work  with 
Dr.  Spring  on  integral  representations  for  the  viscoelastic 
deformation  of  ice.  It  is  concerned  with  different  structures 
of  integral  representation,  and  how  the  respective  kernels 
are  correlated  with  conventional  constant  stress  and  constant 
strain-rate  data. 

Part  II  is  a  paper  written  by  Dr.  Williams,  which  will 
be  published  by  Cold  Regions  Science  and  Technology  during 
1984.  Dr.  Williams  developed  an  alternative  integral 
representation  which  reflected  the  non-monotonic  nature  of 
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the  response  to  constant  stress  without  requiring  an  explicit 
non-monotonic  kernel  function  of  time.  The  latter  is  shown 
to  have  unsatisfactory  physical  implications  as  well  as 
previously  noted  application  difficulties.  The  new 
representation  is  also  shown  to  be  consistent  with  a  variety 
of  different  response  features. 
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Integral  Representations  for  the  Viscoelastic 
Deformation  of  Ice 


U.  Spring  and  L.  W.  Morland 


School  of  Mathematics  and  Physics 
Univesity  of  East  Anglia 
Norwich,  U.K. 


Abstract 

Various  single  integral  representations  which  describe 
non-linear  viscoelastic  response  are  examined  with  regard  to 
the  types  of  test  required  to  determine  the  respective  kernels. 

A  strain  formulation  determined  by  constant  uni-axial  stress 
response  typical  of  ice,  and  its  predictions  for  constant 
strain-rate  response,  are  reviewed,  showing  that  the  latter 
are  sensitive  to  kernel  detail.  An  alternative  stress 
formulation  which  is  determined  by  constant  strain-rate  response 
is  constructed,  and  it  is  shown  that  the  predicted  strain  and 
strain-rate  responses  at  constant  stress  are  compatible  with 
the  typical  responses  exhibited  by  ice. 
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1.  Introduction 

The  viscoelastic  properties  of  ice  are  commonly  exhibited 
by  the  creep  response  at  constant  stress,  and  the  stress 
response  at  constant  strain-rate,  in  uni-axial  compressive 
stress  tests  at  constant  temperature,  reviewed  by  Mellor  (1980). 
While  detailed  results  over  the  large  temperature  range  of 
practical  concern  (230K  -►  273K)  have  not  been  determined. 


the  main  qualitative  features  of  the  responses  up  to  strains 

of  a  few  percent  are  known.  Typical  responses  are  shown  in 

Fig.  1  where  a  is  the  uni-axial  compressive  stress,  e  is 

the  axial  engineering  compressive  strain  defined  as  the  decrease 

in  length  per  unit  initial  length,  t  denotes  time,  and  e 

is  the  strain-rate.  At  constant  stress  a  ,  Fig.  la,  there  is 

where 

an  initial  elastic  strain  jump  ee(o)  =  ct/Eq,  ,  Eq  is  the 

10  -2 

Young's  modulus  at  zero  stress,  of  order  10  Nm  ,  followed 
by  a  primary  decelerating  creep  (e  <  0)  until  time  tm(a) 
when  e  =  0,  e  =  rm(a),  e  =  em(a),  then  a  tertiary  accelerating 
creep  (e  >  0) .  The  minimum  strain-rate  rm(a),  also  referred 
to  as  secondary  or  steady  state  creep,  is  significantly  non-linear 


in  a  and  highly  temperature  dependent,  but  the  strain  em(°) 

6  *“  2 

is  approximately  0*01  over  a  moderate  stress  range  0  -*■  10  Nm 
By  comparison,  the  initial  strain  jump  eg  <<  em  over  this 
stress  range,  and  so  is  negligible  in  applications  involving 
tertiary  creep.  The  associated  strain-rate  shown  in  Fig.  lb 
illustrates  more  clearly  the  non-monotonic  response  features 
not  common  in  other  rheological  models.  At  constant  strain-rate 
e  *  r,  Fig.  lc,  the  stress  increases  continuously  from  zero 
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to  a  maximum  aM(r)  at  time  tM(r)  an<^  strain 
eM(r)  =  rt^r)  ,  then  decreases.  The  strain  eM(r)  is 
also  found  to  be  approximately  0*01  over  a  moderate  range 
of  r  . 

These  multi-type  test  responses  have  been  modelled  by 
viscoelastic  differential  operator  relations  of  fluid  type 
(Morland  and  Spring  1981)  involving  stress,  stress-rate, 
strain-rate,  and  strain-acceleration  terms,  and  of  solid 
type  (Spring  and  Morland  1982)  involving  stress,  stress-rate, 
strain,  and  strain-rate  terms.  In  the  fluid  relation  the 
constant  stress  and  constant  strain-rate  responses  are  neither 
completely  independent  nor  completely  dependent,  so  reflect 
some  common  and  some  distinct  properties.  In  the  solid 
relation  these  two  responses  are  completely  independent,  and, 
in  fact,  are  insufficient  to  determine  the  model  without 
further  reductions.  Mellor  (1980)  conjectured  that  the 
two  types  of  response  should  be  dependent,  which  is  the  familiar 
linear  result.  While  the  differential  fluid  relation  is  closer 
to  this  conjecture,  it  cannot  model  anisotropy  or  strain 
jumps.  In  both  types  it  is  necessary  to  incorporate  dependence 
of  one  or  more  response  coefficients  on  both  stress  and 
strain-rate  or  strain  respectively,  and  covering  adequate  stress- 
strain-rate  or  stress-strain  domains  by  practical  test 
programmes  is  non-trivial  in  uni-axial  stress,  and  the 
difficulty  is  compounded  in  the  multi-axial  tests  required 
to  describe  a  three-dimensional  response. 

In  contrast,  an  integral  representation  expresses  the 
stress  (stress  formulation)  explicitly  in  terms  of  strain 
history  or  the  strain  (strain  formulation)  explicitly  in  terms  of 
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^.history.  The  material  properties  are  therefore  defined  by 
integral  kernel  functions  of  strain  history  or  of  stress 
history,  but  not  both  together.  With  the  above  differential 
relations  (Morland  and  Spring  1981,  Spring  and  Morland  1982), 
the  strain-rate  (strain)  response  to  an  applied  stress 
history,  and  the  stress  response  to  an  applied  strain-rate 
(strain)  history,  are  determined  by  solutions  to  differential 
equations  in  time,  while  with  integral  relations  they  are 
determined  by  direct  integration  or  by  solutions  to  Volterra 
integral  equations.  Integral  representations  may  therefore 
offer  numerical  advantages  in  the  solution  of  boundary-value 
problems,  by  not  increasing  the  order  of  stress  derivatives, 
as  well  as  a  more  direct  correlation  with  test  data.  An 
investigation  of  different  single  integral  representations  to 
describe  given  constant  uni-axial  stress  response  has  been 
presented  by  Morland  and  Spring  (1982) .  Such  strain  formulations 
define  the  uni-axial  response  completely,  and  in  particular, 
the  stress  histories  required  to  maintain  constant  strain-rate 
are  determined  by  Volterra  integral  equations.  Numerical 
solutions  based  on  the  shape  of  respone  illustrated  in 
Figs  la  and  lb  show  that  the  constant  strain-rate  response 
is  acceptable  only  for  one  of  the  representations,  and  is  then 
sensitive  to  kernel  detail,  not  always  giving  the  required 
response  shape  illustrated  in  Fig.  lc.  One  conclusion  is  that 
the  single  integral  strain  representation  should  be  correlated 
approximately  with  both  constant  stress  and  constant  strain-rate 
responses  by  some  optimisation  procedure,  and  not  correlated 
exactly  with  the  constant  stress  response,  but  such  a  scheme 
has  yet  to  be  devised. 


The  above  results  suggest  that  the  constant  strain-rate 
response  of  ice  may  provide  a  better  basis  for  a  single 
integral  representation.  Since  the  kernels  of  the  stress 
and  strain  formulations  discussed  by  Morland  and  Spring  (1982) , 
depending  on  strain  history  and  stress  history  respectively, 
cannot  be  determined  by  constant  strain-rate  response,  we 
now  investigate  two  alternative  stress  formulations  with 
kernels  depending  on  strain-rate  history.  One  defines  a 
viscoelastic  solid  which  exhibits  no  strain  jump  when  a  stress 
jump  is  applied,  which  is  a  good  approximation  in  many 
situations,  and  the  second  defines  a  viscoelastic  fluid.  A 
stress  formulation  is  also  the  most  convenient  for  substitution 
in  momentum  or  equilibrium  equations.  Both  respresentations 
are  equivalent  in  the  small  strain,  small  rotation, 
approximation.  Kernels  describing  the  shape  of  uni-axial 
stress  response  illustrated  in  Fig.  lc  are  constructed, 
and  the  associated  creep  and  strain-rate  responses  at 
constant  stress  are  determined  numerically  from  the  resulting 
Volterra  integral  equation.  It  is  found  that  these  responses 
are  consistent  with  the  shapes  illustrated  in  Figs  la  and  lb, 
and  we  conclude  that  this  form  of  single  integral 
representation  is  a  better  description  of  the  uni-axial 
multi-type  responses  illustrated  in  Fig.  1  than  the  previous 
strain  formulation  determined  by  constant  stress  response. 

2.  Integral  representations 


viscoelastic  material  is 


o(t)  =  RT(t)a  (t)R(t)  =5te<t)] 


(1) 


where  a  is  the  Cauchy  stress,  E  is  a  strain  defined  by 
,T, 


% ( F  F  —  1),  and  F  and  R  are  the  deformation  gradient 
and  rotation  from  a  given  reference  configuration  {see 

A 

for  example,  Dill  (1975)}.  That  is,  the  rotated  stress  a 
at  time  t  depends  on  the  history  of  strain  E(t), 


-  00  <  t  <  t  ,  which  we  designate  a  stress  formulation. 

the  relation 


It  is  supposed  thatjjl)  can  be  inverted  for  viscoelastic 
materials  (thought  not  for  viscous  fluids)  to  give  a  strain 
formulation 

E  ( t)  =  sKito  (x)  ] 


(2) 


in  which  the  strain  at  time  t  depends  on  the  history  of 
the  rotated  stress.  Any  material  symmetries  of  the  reference 
configuration  impose  restrictions  on  the  tensor  functionals 
and  .  If  the  material  is  incompressible,  then  the 
strain  satisfies  the  constraint 


det ( 1  +  2E)  =  1  , 


(3) 


and  the  stress  is  determined  only  within  an  arbitrary 
additive  isotropic  pressure;  that  is,  a  is  replaced  by 

A  , 

o  +  pl  in  the  relations  (1)  and  (2)  where  p  is  an  arbitrary 
scalar  field. 

Multiple  integral  expansions  of  the  stress  functional 


with  kernels  depending  only  on  the  time  lapse  t  -  r 


were  first  proposed  by  Green  and  Rivlin  (1957),  and  have 
received  considerable  attention.  Even  with  the  commonly 
adopted  third  order  truncation  and  restriction  to  isotropic 
material,  the  necessary  test  programme  to  determine  the  kernels 
is  prohibitive.  A  detailed  account  is  given  by  Lockett  (1972). 
While  a  truncated  expansion  is  an  exact  frame  indifferent 
theory,  and  may  describe  highly  non-linear  response  of  a  real 
material,  its  motivation  is  that  the  current  strain  and 
weighted  integration  of  past  strain  (fading  memory)  are 
small.  Furthermore,  the  use  of  multiple  integral  expansions 
in  the  solution  of  boundary-value  problems  is  not  too 
appealing,  particularly  if  a  high  order  truncation  is  required. 
An  alternative  proposal  by  Pipkin  and  Rogers  (1968)  is  that 
the  first  single  integral  term  should  be  a  good  approximation, 
and  reflect  any  significant  non-linearity  by  including  general 
dependence  of  the  kernel  on  the  past  strain  in  a  stress 
formulation  (1),  and  on  the  past  stress  in  a  strain  formulation 
(2) .  Such  single  integral  representations  are 


a ( t)  =  GCE(t) ,0]  + 

A 

E ( t)  =  J[a ( t) ,0]  + 


■  t  t 

GCE(t) , t-x ]dx  , 
J  ~  ~ 

—  00 

'  t  #  /s 

j[a(T) , t-x ]dx  , 

—  00 


(4) 

(5) 


where  the  relaxation  function  G[E,t]  and  creep  function 

/V  A 

J[a,t]  are  tensor  functions  of  strain  E  and  stress  a 


respectively,  and  of  a  positive  scalar  time  argument.  G,  J 
denote  derivatives  of  G  and  J  respectively  with  respect 
to  the  second  (scalar)  argument. 
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It  is  supposed  that  the  stress  after  maintained  zero 
strain  is  zero,  and  the  strain  after  maintained  zero  stress 
is  zero,  so 


% 

.*» 
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G (0, t)  =  0  ,  J(0,t)  =  0  .  (6) 


The  stress  given  by  the  representation  (4)  for  a  constant 
strain  history  E(r)  =  EqH(x)  ,  where  H(t)  is  the 
Heaviside  step  function,  is 


ij 


hi 


E(T)  =  EQH(T)  : 


a  (t)  =  G[Eo,t]  , 


(7) 


and  the  strain  given  by  the  representation  (5)  for  a  constant 

A  A 

(rotated)  stress  history  o(x)  =  oqH(t)  is 


M 


V  , 


a(T)  =  aQH(T) 


E(t)  =  J[oq, t]  . 


(8) 


Thus,  a  programme  of  constant  strain  tests  determines 
the  kernel  G(EQ,t)  over  the  range  of  EQ  covered,  and 
hence  determines  the  single  integral  representation  (4) 

A 

completely.  The  stress  a(t)  then  predicted  for  arbitrary 
strain  history  E(t),  -  00  <  t  <  t  ,  is  given  by  direct 

integration,  and  the  strain  E(t)  predicted  for  arbitrary 

A 

stress  history  o(t)  ,  -  <®  <  x  <  t  ,  is  the  solution  of  a 
non-linear  Volterra  integral  equation.  Similarly,  a  programme 

A 

of  constant  stress  tests  determine  the  kernel  J(aQ,t)  over 
the  range  of  aQ  covered,  and  hence  determines  the  single 
integral  representation  (5)  completely.  Pipkin  and  Rogers 
(1968'  show  how  a  sequence  of  multiple  integral  corrections 
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can  be  constructed  if  the  responses  to  multi-step  strain  or 
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stress  tests  do  not  correspond  exactly  to  the  predictions  of 
the  lower  order  representations,  starting  with  the  relations 
(4)  and  (5)  respectively,  but  this  is  again  a  prohibitive 
programme  (Lockett  1972).  Clearly  single  integral  representations 
are  the  most  tractable  for  stress  and  deformation  analysis,  and 
Pipkin  and  Rogers  (1968)  have  demonstrated  for  some  materials 
that  the  strain  formulation  (5)  in  uni-axial  stress  determined 
by  a  single-step  constant  stress  test  predicts  good  agreement 
with  the  responses  to  multi-step  loads.  If,  however, 
multi-type  test  data  are  more  common,  such  as  the  constant 
stress  and  constant  strain-rate  responses  of  ice  illustrated 
in  Fig.  1,  there  is  no  direct  scheme  analogous  to  the 
multi-step  test  case  to  correct  a  single  integral  approximation. 
Morland  and  Spring  (1982)  have  investigated  the  strain 
formulation  (5)  in  uni-axial  stress  with  varying  scalar 
kernels  J[a,t]  compatible  with  the  shape  of  the  small  strain 
creep  responses  Figs  la  and  lb.  Numerical  solutions  of  the 
non-linear  Volterra  integral  equations  for  the  corresponding 
stress  histories  a(t)  required  to  maintain  constant 
strain-rate  in  t  >  0  show  that  the  required  shape  of 
response  Fig.  lc  is  not  always  obtained.  That  is,  the  response 
to  constant  strain-rate  predicted  by  a  single  integral 
representation  (5)  is  sensitive  to  kernel  variation,  which 
implies  that  a  kernel  J[ o , t]  determined  solely  by  constant 
stress  data  may  yield  unsatisfactry  predictions  for  other 
loading  responses.  The  representation  (5)  should  therefore 
be  correlated  approximately  with  both  constant  stress  and 
constant  strain-rate  responses  in  some  optimal  manner,  but  no 
direct  scheme  has  yet  been  devised. 
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The  above  results  suggest  that  the  constant  strain-rate 
responses  of  ice  may  provide  a  better  basis  for  a  single 
integral  representation.  In  the  stress  formulation  (4), 
constant  strain-rate  implies  a  varying  argument  E(t)  in 
the  kernel  G  [E(T),t-r]  ,  and  no  unique  correlation  with 
the  function  G[E,t];  only  constant  strain  response  determines 
G[E,t].  Alternately/  in  the  strain  formulation  (5),  a 

A 

non-uniform  response  a(t)  to  constant  strain-rate  in  the 

A 

kernel  j[a (t) ,t-i]  yields  a  sequence  of  integral  equations 
for  J  with  no  unique  solution.  The  same  situation  arises 
if  constant  strain  is  prescribed  in  the  representation  (5) , 
discussed  in  more  detail  by  Morland  and  Spring  (1982) ,  which 
shows  that  G  does  not  determine  J  ,  and  similarly  J  does 
not  determine  G  .  That  is,  the  single  integral  stress  and 
strain  formulations  (4)  and  (5)  are  not  equivalent,  in  contrast 
to  the  relaxation  and  creep  integral  relations  of  linear 
viscoelasticity  for  which  the  relaxation  and  creep  functions 
(of  time  only)  are  related  uniquely  by  a  linear  Volterra 
integral  equation.  Hence,  a  kernel  which  can  be  determined 
by  constant  strain-rate  response  must  depend  only  on 
strain-rate  history  and  time,  so  we  now  examine  two  such 
integral  representations  of  the  stress  formulation  (1) . 

The  direct  analogy  to  the  representation  (4)  is 
-  ft  . 

a (t)  =  T  C  E '  ( t)  ,0]  +  T  CE'  (t)  ,t-T]dx  ,  (9) 

—  oo 

with 


T  (0,t) 


0 


(10) 


so  that  the  response  to  constant  strain-rate  E'(t)  =  rsH  (x  ) 


E'  (t)  =  rsH(x)  : 


o  (t)  =  TsCrs, t]  , 


which  determines  the  kernel  T  [r  ,  tl  over  the  range  of  r 

~  s  -  s  ~  s 


covered.  Since  E'(t)  is  a  rate  of  strain  relative  to  a 


given  reference  configuration,  the  representation  (9)  defines 


a  viscoelastic  solid,  and  since  bounded  argument  E'(t)  is 


implied,  no  strain  jumps  arise;  a  stress  jump  induces  a  jump 


in  strain-rate.  This  is  a  satisfactory  approximation  for  ice 


in  many  situations.  Also,  anisotropy  of  the  reference 


configuration  can  be  described  through  the  structure  of  the 


tensor  kernel  T  [E',t], 


An  alternative  viscoelastic  fluid  relation  can  be  constructed 


from  Noll's  reduced  form  of  the  representation  (1)  when 


dependence  on  the  strain  E  and  rotation  R  from  a  given 


reference  configuration  is  eliminated,  namely 


o(t)  =  fiCEt(x)]  , 


where  is  an  isotropic  tensor  function  of  the  tensor  Efc  to 


satisfy  frame  indifference.  Here  Efc(x)  is  the  strain  at  time 


x  relative  to  the  configuration  at  time  t  ,  and  is  independent 


of  any  given  reference  configuration.  All  configurations  are 


necessari 


ly  isotropic.  Et'(x)  is  the  strain-rate  at  time  x 


relative  to  the  current  configuration.  The  appropriate  integral  !*] 


representation  is 


Tf[Ft'  (T) ,t-T]dT  , 


(13) 


o(t)  =  TfCEt'  (t)  ,0]  +  j 

—  00 


which  satisfies  the  frame  indifferent  requirement  (isotropy) 
if  Tf  has  a  Rivlin-Ericksen  expansion 

Tf  [  Ef '  (t ) , t]  =  <j»Ql  +  cJ,1Et  *  (t  )  +  4>2{Et'(T)}2  (14) 


-■■4 


•M 


where  the  response  coefficients  4>a  (a  =  0,1,2)  are  functions 
of  the  invariants  of  E^' (t)  and  time.  Also 

Tf(0,t)  =  0  <=>  4>q  =  0  when  Et'(T)  =  0.  (15) 


Thus  Tf  is  a  polynominal  in  the  tensor  Et'(x),  but  has 
arbitrary  dependence  on  its  invariants.  If  the  solid  defined 
by  the  representation  (9)  is  isotropic  in  the  reference 
configuration,  then  E  is  replaced  by  E  =  H (FF  -  1)  , 

A 

o  by  o  ,  and  Tg  has  an  expansion  in  E'  analogous  to  the 
expansion  (14)  of  Tf  in  Et’’  The  response  to  constant 
strain-rate  Et'(i)  =  rfH(x)  given  by  the  representation  (13) 
is 

Et'(i)  =  rf  H  (  t)  :  a(t)  =  Tf[rf,t]  (16) 


which  determines  the  kernel  Tf[rf,t]  over  the  range  of  rf 
covered. 

While  the  strain-rate  D(t)  relative  to  the  current 
configuration,  defined  by  the  symmetric  part  of  the  spatial 
velocity  gradient  at  time  t  ,  is  simply 

D ( t)  =  E.  '  (t)  ,  (17) 


a  kernel  dependence  on  D(t)  would  require  rates  measured  with 
respect  to  the  continuously  changing  sequence  of  configurations 
as  the  integration  variable  t  changes,  and  does  not  appear 
to  be  a  natural,  or  practical,  representation.  Also  the  fluid 
representation  (13)  ,  (14)  is  distinct  from  the  isotropic  form 

of  the  solid  relation  (9) ,  since  the  relative  strain-rate 
Et'(x)  and  strain-rate  E'(x)  from  a  reference  configuration 
are  not  identical.  However,  for  small  strain  | |e| |  <<  1  and 
small  rotation  | |R  -  l| |  <<  1  ,  if  terms  of  these  magnitudes 
are  neglected  in  comparison  with  unity,  then 

E*  (x)  =  E'  (x)  =  E t '  (x)  .  (18) 

In  this  situation  the  constant  strain-rates  r  and  r c  in 

~s 

the  results  (11)  and  (16)  are  the  same,  and  the  kernels  T 

~  b 

and  Tf  are  the  same,  so  the  representations  (9)  and  (13)  are 
equivalent.  If  though  a  common  function  T  is  adopted  for  Tg 
and  Tf  and  extended  to  large  strain,  with  small  or  large 
rotation,  the  respective  arguments  E'(x)  and  Et'(x)  become 
different  and  the  representations  (9)  and  (13)  are  distinct. 

We  will  adopt  the  common  representation  for  small  strain  and 
zero  rotation  to  investigate  the  uni-axial  stress  description 
determined  by  a  typical  response  illustrated  in  Fig.  lc . 

3.  Uni-axial  stress  relations 

For  uni-axial  compressive  stress  o  ,  axial  compressive 
engineering  strain  e  and  equal  lateral  compressive  engineering 
strains  e^,  zero  rotation,  and  the  incompressibility 


approximation, 


Gn  =  “a  '  G^  -  0  (*■¥  i.  1  *  i). 


ij 


(19) 


E11  =  -e(l  -  *se)  ,  E22  =  E33  =  *je(l  -  e) 


-1 


In  the  small  strain  approximation  E^  =  -  e  ,  both 
representations  (9)  and  (13)  yield  an  axial  relation 


a (t)  =  T[e ' (t) ,0]  + 


T[e '  (t ) , t-r ]dr 


(20) 


where  T[r,t]  is  the  stress  a(t)  required  to  maintain  a 
constant  strain-rate  eft)  =  rH(t).  Later  calculations  will 
use  functions  T[r,t]  constructed  from  families  of  stress 
responses  for  different  r  with  the  shape  illustrated  in 
Fig.  lc.  Given  the  kernel  T(r,t)  and  a  prescribed  stress 
history  a(t),  the  relation  (20)  is  a  non-linear  Volterra 
integral  equation  for  the  strain-response  e(t).  Incompressibility 
implies  that  the  fluid  relation  (13)  and  isotropic  form  of  the 
solid  relation  (9)  have  only  two  independent  principal  deviatoric 
components,  since  the  deviatoric  kernel  T  satisfies  tr  T =  0  , 
but  the  uni-axial  configuration  (19)  yields  only  one 
independent  component  with 


T22  “  2  T11 


(21) 


Thus  tests  with  two  independent  stress  components  are  required 
to  determine  the  deviatoric  tensor  relation. 

Let  a  constant  stress  a(t)  =  acH(t)  be  prescribed  in 
the  integral  equation  (20),  so  e(t)  =  0  for  t  £  0. 


Define 


fci  =  i6  ,  ri  =  e'(ti),  i  =  0,1,2,...  ,  (22) 

where  6  is  a  small  time  interval,  and  interpret  rQ  as  the 
strain-rate  at  t  =  0+;  that  is,  the  strain-rate  jump  at 
t  =  0.  By  equation  (20), 

•  a  c  =  T[rQ,0],  (23) 

which  has  a  unique  solution  rQ(ac)  if  T[r,0]  is  monotonic 
increasing  in  r;  that  is,  if  constant  strain-rate  rH(t) 
requires  an  initial  stress  jump  a  which  increases  with  r  , 
assumed  here.  This  stress  jump,  instead  of  the  smooth  rapid 
rise  from  zero  illustrated  in  Fig.  lc,  is  a  necessary  feature 
of  representations  with  no  strain  dependence,  and  will  be 
incorporated  in  the  model  functions.  Now  apply  the  trapezoidal 
rule  to  the  first  time  interval  of  equation  (20),  then 

°c  =  T[ri'0]  +  %fi{TCrL,0]  +  T[  rQ,  6  ]  }  (24) 

which  is  an  implicit  equation  for  r^.  Again,  assuming  that 
T[r,0]  is  monotonic  increasing  in  r  ;  that  is,  the  initial 
stress-rate  at  constant  strain-rate  r  increases  with  r  ,  there 
is  a  unique  solution  (a  ) .  Continued  application  of  the 
trapezoidal  rule  over  successive  time  intervals  gives  for 
n  ~~  2,3,... 

°c  =  TCrn'0]  +  %fiT[rn,0] 

.  n-1  . 

+  5{»sT[r  ,t  ]  +  Z  T[r.  ,t  ,]}, 
on  •;— i  iiii. 


(25) 


once 


which  is  an  implicit  equation  for  rn 

(i  =  0,1,..., n-1)  are  determined,  with  unique  solution 
by  the  above  monotonicity  properties. 

Although  the  algorithm  (23)  -  (25)  yields  a  sequence 
of  implicit  equations  for  each  rn  (n  =  1,2,...),  the 
unknown  rn  occurs  only  in  a  monotonic  function  comprising 
two  terms,  and  numerical  solution  is  simple  and  accurate. 

Tests  of  accuracy,  and  the  main  solution  features,  can  be  made 
by  comparison  with  exact  solutions  for  an  associated  linear 
viscoelastic  relation  given  by  taking 

T[e '  ,  t]  =  T^tje'  ,  TCr^t^]  =  TJl(tj)ri,  (26) 

when  the  axial  relation  (20)  becomes 

a(t)  =TA(0)e'(t)  +|  (t-r)e’ (T)dr  .  (27) 

— oo 

Let  the  Laplace  transform  of  f (t)  defined  on  t  ^  0  be 
denoted  by  L[f(t)]  =  f(s),  and  define  f(s)  =  sf(s),  then 
for  strain  and  stress  histories  zero  on  t  <  0  ,  the  transform 
of  the  convolution  relation  (27)  gives 

e7 (s)  =a(s)/T^(s).  (28) 

For  the  illustrations,  T^(t)  is  represented  by  a  series  of 
exponentials,  and  for  a(s)  =  o^/s  ,  e'(s)  can  then  be  inverted 

analytically  to  express  e'(t)  as  a  series  of  exponentials. 

4.  Illustrations 

The  constant  strain-rate  response  illustrated  in  Fig.  lc, 
but  with  an  initial  non-zero  stress,  has  been  modelled  by  the 


uiiu 


„•  -■  ■>  V 


-  17  - 


kernel  functions 


T[r,t]  =  cM(r){A  +  Be  ^rt  -  De  drt} 


for  three  different  sets  of  constants  A,B,D,b,d,  with 

d  >  b  >  0.  Since  t  enters  only  in  the  product  e  =  rt  ,  the 
maximum  stress  for  each  r  occurs  at  a  strain  eM  given  by 


exp{ (d-b)eM>  =  dD/bB  , 


which  is  indpendent  of  r  ,  consistent  with  observations.  We 
adopt  the  value  eM  =  0*01,  then  for  maximum  stress  • 


A  .  D  -0*01b  _  -0*0ld  . 

A  +  Be  -  De  =  1  . 


The  initial  stress  is 


T[r,0j  =  oM(r){A  +  B  -  D}  , 


which  requires 


A  +  B  -  D  >  0  , 


and  a  bounded  positive  long  time  (large  strain)  limit 
AoM(r)  <  oM(r)  requires 


1  >  A  >  0  . 


Also,  the  strain  e^.  (>  eM)  at  inflexion,  T  =  0,  is  given  by 


exp{(d-b)eI)  =  d2D/b2B 


(35) 


We  prescribe  A,  d/b,  Sj/eM  with  e^  =  0*01  to 
calculate  the  three  sets  of  constants  shown  in  Table  1  from 
the  relations  (30),  (31),  (35). 


I 

II 

III 

eI/eM 

1*75 

1-84 

1-90 

A 

0*7 

0-7 

0-7 

B 

0-6412 

0-4520 

0-3191 

D 

1-0965 

0-7009 

0-5323 

b 

53-6479 

30-4575 

5-1685 

d/b 

5 

10 

100 

Table  1.  Parameters  for  three  families  of  constant 
strain-rate  responses 

Figure  2  shows  the  corresponding  curves  for  T[r,t]/a  (r)  =  a/a 

“  M 

as  functions  of  e  =  rt.  Following  Mellor's  (1980)  conjecture, 

the  maximum  stress  relation  c,.(r)  is  assumed  to  be  the 

inverse  function  of  the  minimum  strain-rate  rm(a)  at  constant 

stress,  and  we  adopt  a  close  least  squares  fit  to  laboratory 

data  at  273  »13K  over  the  axial  stress  range 
5  -2 

0  s  04  9  *  10  Nm  derived  by  Smith  and  Morland  (1981): 

oM(r)  =  1 • 991 tan~ 1 ( 1 • 609r )  +  4 •2929tan~1 (0-0206r) ,  (36) 

5-2  -1 

where  a  is  measured  in  units  10  Nm  and  r  in  units  a 
o^' (r)  >0  so  T[r,0]  is  monotonic  increasing  as  required, 
and  further,  T[r,0]  is  monotonic  increasing  as  required  since 


dD  >  bB  by  equation  (30)  with  d  >  b  . 


by 


A  comparison  linear  viscoelastic  relation  (27)  is  given 


T£(t)  =  k{  A  +  Be”bt  -  De_dt} , 


(37) 


where  k  is  a  constant.  Figure  3  shows  the  stress 

o  =  T^(t)r  as  a  function  of  e  =  rt  for  three  values  of 

constant  strain-rate  r  using  the  parameter  set  I  in  Table  1 

5  -2-1 

and  k  =  4*5  x  10  Nm  a  .  Here  the  maximum  stress  occurs 

at  constant  t„  =  0*0la  for  all  r  ,  and  hence  at  e.,  =  rt., 
M  MM 

which  is  proportional  to  r.  The  Laplace  transform  solution 
(28)  for  o ( t )  =  a  H(t)  is 

a  (s+b) (s+d) 

e'  (s)  =  - 2 - 

ks (A (s+b) (s+d)  +  Bs(s+d)  -  Ds(s+b)} 


°_c  U _ V_  +  _W_  1 

k  (s  s+v  s+w  j  ' 


(38) 


giving 


e' (t)  =  ~  H(t)|u  -  Ve~Vt  +  We"wt j 


(39) 


With  the  above  parameters  and  t  in  units  a 


U  =  1*4286,  V  =  0*5608,  W  =  3*2183,  v  =  30*428,  w  =  1352*71, 


and  Fig.  4  shows  the  corresponding  ke'(t)/ac  .  Comparison 
values  at  different  times,  including  the  time  ^m^0c^  =  36>*(>hrs 
to  minimum  strain-rate,  of  the  analytic  solution  ke'(t)/ac  ant 
that  given  by  the  numerical  algorithm  (23)  -  (25)  applied  to 


the  integral  equation  (27),  are  shown  in  Table  2.  The  numerical 


t  hrs 

ke ' ( t) /o  (anal ) 

ke ' ( t) /a  (num) 

0 

4*0861 

4*0861 

10 

1*5740 

1*5739 

36-6 

0*9460 

0*9460 

100 

1*0323 

1*0323 

1000 

1*4112 

1*4105 

error 

0 

-1  x  lo' 

O 

O 

-7  x  io' 


Table  2.  Comparison  of  analytic  and  numerical  strain-rates 

at  constant  stress  for  linear  viscoelastic  relation. 


algorithm  is  clearly  accurate. 

Applying  the  algorithm  to  the  non-linear  relation  (29) 

with  parameter  sets  I,  II,  III  in  turn  gives  'he  respective 

strain-rates  at  different  constant  stresses  a  .  The 

c 

minimum  strain-rate  r  (o  ) ,  time  to  minimum  t  (a  ) ,  and  strain 

me  me 

at  minimum  e  (o  )  ,  are  given  in  Table  3.  The  strain  e  (a  ) 
me  me 
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0  10  Nm 
c 

,  ^  “I 

r  (o  )  a 
m  c 

tm(0c,hrs 

e  (a) 
m  c 

0*25 

0*0806 

1048 

0*0111 

0*5 

0*1542 

524 

0*0114 

T 

0*75 

0*2540 

348 

0*0119 

X 

1*0 

0*3549 

258 

0*0128 

1*5 

0*6193 

165 

0*0155 

2*0 

1*0863 

109 

0*0196 

II 

0*5 

0*1597 

524 

0*0107 

1*0 

0*3429 

254 

0*0113 

III 

0*5 

0*1560 

628 

0*0121 

1*0 

0*3326 

306 

0*0127 

Table  3.  Constant  stress  results  for  three  Darameter  sets. 


is  approximately  0*01  for  all  three  sets  for  oc  <  10  Nm  , 

5  -2 

but  increases  to  approximately  0*02  at  ac  =  2  *  10  Nm 

Figure  5  shows  the  normalised  strain-rate  e'(t)/rm(ac) 

as  a  function  of  t/t  (o  )  for  different  values  of  o  . 

me  c 

5  -  2 

The  response  at  ac  =  2  *  10  Nm  is  not  satisfactory  for  long 

times,  but  the  corresponding  creep  curves  showing 
•t 

e(t)  =  e' (t)dr  in  Fig.  6  indicate  that  the  strain  e  is 
o 

becoming  large  and  outside  the  scope  of  the  representation.  The 

responses  at  low  stresses  are  very  similar,  and  are  consistent 

with  the  required  shapes  in  Figs  lb  and  la.  Figures  7  and  8 

5  -2 

compare  the  strain-rate  and  creep  responses  at  oc  =  10  Nm 
for  the  three  parameter  sets,  showing  that  the  response  is  not 
sensitive  to  kernel  detail  when  the  representation  and  kernel 
are  determined  by  constant  strain-rate. 

5.  Conclusions 

Single  integral  representations  of  non-linear  viscoelastic 
response  are  determined  fully  by  one  type  of  test,  but  different 
forms  of  representation  are  required  for  different  types  of 
test,  and  they  are  not  equivalent.  The  integral  kernel  is 
determined  directly  by  the  appropriate  test  response,  in  one, 
two,  or  three-dimensional  tests,  and  provides  a  much  simpler 
correlation  with  test  data  than  the  viscoelastic  relations 
of  differential  type.  An  integral  representation  also  offers 
numerical  advantages  in  the  solution  of  boundary-value  problems. 
It  has  been  shown  that  stress  formulations  for  ice  depending 
only  on  strain-rate  history,  not  strain,  of  solid  type  and  of 


fluid  type,  which  are  determined  by  the  stress  history  required 
to  maintain  constant  strain-rate,  can  predict  an  appropriate 
shape  of  strain-rate  response  to  constant  uni-axial  stress 
for  small  strain.  A  strain  formulation  based  on  constant 
stress  response  predicts  less  satisfactory  response  to  constant 
strain-rate.  An  improved  single  integral  approximation  may 
be  possible  by  correlating  the  above  stress  formulation 
approximately  with  constant  strain-rate  data  within  a  range 
of  kernel  parameters,  and  seeking  an  optimum  parameter  set 
for  approximate  agreement  with  constant  stress  data. 
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Figure  Captions 


1.  Typical  responses  of  ice  in  uni-axial  compressive 
stress  a  : 

(a)  strain  e(t)  at  constant  stress,  (b)  strain-rate 
e(t)  at  constant  stress ,  (c)  stress  a (t)  at  constant 
strain-rate . 

2.  Constant  strain-rate  responses  for  parameter  sets 

1  ( - )  ,  II  ( - )  ,  and  III  ( - )  . 

3.  Linear  viscoelastic  stress  response  at  constant  strain-rate 

4.  Linear  viscoelastic  strain-rate  at  constant  stress. 

5.  Strain-rate  response  to  constant  stresses  for  parameter 

set  I:  a  =  0*25  (  —  -••),  1  (- - ),  1*5  ( - ), 

2  (- • - • )  x  105Nm-2. 

6.  Strain  response  to  constant  stresses  for  parameter  set 

I:  oc  =  0*25  1( - ),  1*5  ( - ),  2  ( - )  *  l05Nm" 

5  -2 

7.  Strain-rate  response  to  constant  stress  oc  =  10  Nm 

for  parameter  sets  I  ( - ),  II  ( - ),  and  III  (-•-•). 

,  5  -2 

Strain  response  to  constant  stress  oc  =  10  Nm 

parameter  sets  I  ( - ),  II  ( - ),  and  III  (-•-•) 


8. 


for 


a-const 


(c; 


r -const 


e 


Figure  4.  Linear  viscoelastic  strain-rate  at  constant  stress 
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SINGLE  INTEGRAL  REPRESENTATIONS  IN  ICE  MECHANICS 


H.  T.  Williams 

School  of  Mathematics  and  Physics, 
University  of  East  Anglia,  Norwich,  U.K. 


Abstract 

A  single  integral  viscoelastic  constitutive  equation  for 
ice  is  developed  which  possesses  significant  theoretical  and 
practical  advantages  over  previously  suggested  equations  of 
this  type  (Spring  and  Morland,  1983).  The  theory  is  specialised 
to  the  case  of  small  strain  uniaxial  compression  and  the 
resulting  constitutive  equation  is  shown  to  verify  the  relations 
between  the  experimental  data  obtained  in  constant  load  (CL) 
creep  tests  and  constant  displacement  rate  (CD)  "strength" 
tests  conjectured  in  Mellor  and  Cole  (1982)  and  demonstrated 
in  Mellor  and  Cole  (1983) . 
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type,  and  presupposes  that  the  material  under  consideration  is 
mechanically  incompressible  (a  reasonable  approximation  for 
ice) . 

The  three  dimensional  small  strain  representation  thus 
obtained  is  then  further  specialised  to  the  case  of  uniaxial 
compression,  and  three  types  of  uniaxial  compression  test  are 
then  considered.  These  tests  are  the  constant  strain  (stress 
relaxation)  test,  the  constant  strain-rate  test  and  the  constant 
stress  (creep)  test.  The  former  test  is  not  significant 
experimentally  but  is  theoretically  important  because  it  imposes 
physical  conditions  on  the  model  which  must  be  met. 

Exact  solutions  are  obtained  for  the  constant  strain  and 
constant  strain-rate  tests  but  the  solution  of  the  constant 
stress  test  involves  a  non-linear  Volterra  integral  equation  of 
the  first  kind,  which  presents  serious  difficulty  both  numerically 
and  analytically.  This  difficulty  can,  however,  be  circumvented 
by  means  of  an  approximate  solution,  which  bears  out  the 
conjectures  of  Mellor  and  Cole  (1982,  Table  1,  p.  206)  on  the 
correspondence  between  constant  strain-rate  and  constant  stress 
data. 

To  conclude,  a  possible  empirical  form  of  the  kernel  of 
the  single  integral  law  is  examined. 


The  Constitutive  Equation 


As  a  starting  point  consider  the  reformulation 

a  (t)  =  [  J  (t)  ] 

—  00 


(1) 


of  the  constitutive  equation  given  for  the  Cauchy  stress  in 
a  viscoelastic  fluid  by  Spring  and  Morland  (1983,  eq.  (12)), 
where 


Jt(T)  =  F(t)F_1(T) {F-1(t) }T  FT(t) 


(2) 


.(Lodge,  1974,  p.105) 

is  the  Finger  strain  tensor4and  F  is  the  deformation  gradient 
tensor.  Following  Morland  and  Spring  (1982),  suppose  that  the 
constitutive  equation  (1)  can  be  expressed  as  a  single  integral 
representation 


o  (t) 


—  00 


f[Jt(T)  , 


t 


i]dx 


(3) 


While  the  single  integral  formulation  (3)  appears  different  to 
that  adopted  by  Morland  and  Spring  (1982),  it  can  be  shown  that 
the  two  formulations  are  equivalent. 

In  order  to  make  further  progress  it  is  necessary  to  make 
simplifying  assumptions.  Suppose  that  the  material  is 
mechanically  incompressible,  a  conventional  and  satisfactory 


approximation  which  leads  to  the  constraints 


det  J  (t)  =  1  ,  det  F  =  1 


(4) 


and  requires  that  the  constitutive  equation  (3)  only  determines 
the  Cauchy  stress  up  to  an  arbitrary  isotropic  pressure  p  , 
where 


o(t)  =  -  pi  +  oE(t) 


(5) 


oE  is  known  as  the  extra  stress.  It  should  be  noted  that  the 
decomposition  (5)  is  not  unique  but  is  commonly  made  so  by 
taking  c?E  to  be  the  deviatoric  stress:  this  is  not  adopted 
here,  however,  because  it  leads  to  an  inconvenient  form  for 
the  reduced  constitutive  equation  derived  later. 

If  the  constitutive  equation  (3)  is  to  provide  a  satisfactory 
basis  for  a  physical  theory  of  viscoelastic  materials  it  must 
comply  with  the  principle  of  material  frame  indifference 
(objectivity) .  This  requirement,  together  with  the  restrictions 
(4)  and  (5)  imposed  by  incompressibility,  leads  to  the 
Rivlin-Ericksen  representation 

o (t)  =  -  pi  +  [  f, (I  ,11  ,  t  -  t)J  (T)dT 
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IB  =  tr  B  and  IIB  =  tr  B 


(12) 


It  should  be  noted  that,  in  equation  (9) ,  a  term  given  by 

t  t 

|  f x  (3 ,  3,  t  -  x )  I  dx  +  |  f_x(3,  3,  t  -  x )  I  dr 

o  o 


has  been  incorporated  into  the  arbitrary  pressure  p  . 

The  comparison  of  equation  (9)  with  the  well-known 
constitutive  equation 


a  =  -  pi  + 


yl^IB,IIB^? 


V-L  (IB'-' 


IIB)! 


-1 


(13) 


of  an  incompressible  isotropic  elastic  material  provides 
motivation  for  further  model  restrictions.  A  particularly 
simple  version  of  equation  (13)  is  the  Neo-Hookean  constitutive 
equation  obtained  by  setting  =  0  and  taking  to  be 

a  constant.  The  corresponding  version  of  equation  (9)  is 
obtained  in  a  similar  manner  with  the  exception  that  ^  must 
now  be  taken  to  be  a  function  of  t  alone.  It  follows  that, 
for  the  "Neo-Hookean"  viscoelastic  fluid,  the  stress  is  given 
by 


a  (t) 


(t 


-  T  ) Jfc  (T)dT 


(14) 


It  will  emerge  in  section  3  that  (14)  is  incapable  of  providing 
a  satisfactory  explanation  of  the  non-monotonic  response  of  ice 


'  «r  ;  -4~ 
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in  a  constant  strain-rate  test  because  of  its  linear  dependence 
of  stress  on  relative  strain  history.  This  difficulty  is 
avoided  by  the  more  general  version 


a ( t)  =  -pi 


-ir 


K,2CT<IJ 


-  3) 


-  T]Jt(x)dx 


(15) 


where  the  kernel  K  is  defined  by 


Ktjdj  -  3)  ,  t]  =  3y1(IJ,  t) 


(16) 


and  K/;l  and  K,2  denote  the  partial  derivatives  of  the 
kernel  K  with  respect  to  its  first  and  second  arguments, 
respectively. 


3.  Small  Strain  Uniaxial  Tests 

In  this  section  a  one-dimensional  constitutive  equation 
is  derived  from  equation  (15)  which  can  be  applied  to  the 
uniaxial  experimental  tests  discussed  previously.  In  a  uni-axial 
stress  test 


(17) 


The  corresponding  deformation  gradient  history  is 


V  v-  V  VA^  MWTV  VJ'-’V  V-V«  ,v.v  ,V .*  .  -  .’•  .  •  ^  - 


•  _  ■  _  «  _  '  _  " 


F  { t)  =  1  +  e  ( t) 


1  "  *5£(t) 


1  -  ^e(t) 


which  approximately  satisfies  the  incompressibility  constraint 
(4) 2  provided  that  the  small  strain  condition 


e  ( t)  <<  1 


holds  for  the  longitudinal  strain  e  (t) .  The  substitution 
of  (18)  into  the  constitutive  equation  (15)  via  the  relation 
(2)  leads  to  the  two  scalar  equations 


=  "P  -  jj  K,2Cet(x)Z'  "  T]C1  " 


2et (t) ]dx 


°22  °33 


=  -P  "  yj  K, 2Cet(T)2'  ”  T]C1  +  et(T)]dT'  (2: 


whe  re 


efc(x)  =  e  ( x )  -  e  (t) 


The  result 


a(t)  =  [  K-Ce.(x)2,  t  -  xje^fxjdx 
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now  follows  through  the  elimination  of  the  pressure  between 
(17),  (20)  and  (21). 

As  an  example  of  how  the  uniaxial  constitutive  equation 
(23)  can  be  specialised  to  a  specific  test  consider  the 
constant  strain  test  given  by  the  linear  strain  history 


e  (t)  =  0  ,  t  <  0  , 


e  (t)  =  e  ,  t  >  0  , 


The  substitution  of  (24)  into  (22)  and  (23)  gives  the 
expression 


a  (t)  =  K(e“2,t)e  ,  t  >  0  , 


for  the  uniaxial  stress,  assuming  the  physically  expected, 
monotonic,  long-time  behaviour 


K (X,  t)  -*■  0  and  K  P(X,  t)  -*■  0  as  t  -*■  00 


which  must  hold. 


4 .  Constant  Strain-Rate  Test 

The  constant  strain-rate  test  is  specified  by  the  linear 


KO 


strain  history 


e(t)  =0  ,  t  <  0 


e ( t)  =  rt 


V- 


where  r  is  the  strain  rate.  Using  the  alternative  form 


a  ( t) 


rt 

=  (KCe 


(t)2,  t-i]  +  2e  t  (  t) 2K  1[et(x) 2 , t-x]}e^(x)dx  (28) 


of  equation  (23),  obtained  by  integrating  by  parts,  with  the 
dash  denoting  partial  differentiation  with  respect  to  x  , 
the  history  (27)  gives  the  relation 


a(t)  =  j  {K[r2  ( t—  t  )  2 ,  t-x]  +  2rz  (t-x^K^Cr2  ( t-x2  )  ,  t-x]  }rdx  .  (29 
o 


The  results  of  the  constant  strain-rate  test  are  usually 
recorded  as  a  one  parameter  family  of  curves  of  stress  measured 
against  strain  (rt)  ,  the  parameter  being  the  strain  rate. 

Thus  the  experimental  data  is  recorded  as  a  function  of  two 
variables  which  can  be  expressed  in  terms  of  the  model  by 
substituting  u  =  r(t-x)  in  (29)  to  obtain 


^  —  1  —  1 

a(e,  r)  =  [K(u2,  r  u)  +  2uzK,.(u2,  r  u)  ]  du  . 

Jn  1 


The  relation  (30)  can  be  inverted  to  give 


K  (e  2 ,  t)  =  c 


1|  5  .(u,  t-1u 

Jn  '  ± 


)  du  , 


where  a  ,  and  5  -  are,  respectively,  the  partial 

f  ^ 

derivatives  of  a  with  respect  to  its  first  and  second  arguments. 


S  S 


•  ■  '> 
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Thus  the  kernel  can,  in  principle,  be  determined  from  the 
constant  strain- rate  test  by  means  of  equation  (31) . 

When  the  results  given  above  are  specialised  to  the  simple 
model  (14),  the  kernel  is  a  function  of  t  only  and  equation 
(31)  can  be  replaced  by 


K  (t)  =5  .  (rt ,  r)  . 
/ 


It  can  be  shown  that  this  model  is  equivalent  to  the  linear 

viscoelastic  relation  employed  by  Spring  and  Morland  (1983) 

to  test  the  accuracy  of  their  numerical  algorithm  and,  as 

can  be  seen  from  Spring  and  Morland  (1983,  Fig.  3),  the  model 

cannot  adequately  represent  the  behaviour  of  ice  in  this 

test  because  the  strain  at  maximum  stress  is  proportional  to  the 
constant  strain-rate  r  and  the  maximum  stress  is  also 
proportional  to  r,  which 

contradicts  the  experimental  evidence  (see,  for  example, 

Mellor  and  Cole,  1983)  that  the  strain  at  maximum  stress 
(typically  around  1%)  is  largely  independent  of  r  and  the 
maximum  stress  is  a  non-linear  function  of  r  ,  tending  to 
a  finite  value  as  r  -*■  00  .  There  is  a  further,  more  important, 
reason  for  the  failure  of  the  "neo-Hookean"  fluid  model  (14) . 
When  equation  (32)  is  used  to  determine  the  kernel  from  a 
typical  constant  strain  rate  response  (e.g.  Spring  and  Morland, 
1983,  Fig.  3(c)),  K(t)  is  found  to  be  non-monotonic  and 
decreases  through  zero  attaining  a  negative  minimum  value  before 
decaying  to  zero.  This  kernel  is  physically  unacceptable 
because  it  predicts  absurd  reponse  in  the  relaxation  test 
equation  (25)  and  it  is  not  a  monotonic  decreasing  function 
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of  time  (a  reasonable  requirement  for  a  one-dimensional 
relaxation  function,  although  it  is  a  matter  of  some  controversy 
whether  non-monotonic  relaxation  functions  are  compatible 
with  the  principles  of  thermodynamics) . 

Equation  (25)  shows  how  the  present  model  overcomes  the 
difficulties  encountered  above.  Representing  the  constant 
strain  test  data  (if  obtainable)  by  a  function  o*(e,  t)  , 
where 


o* (e ,  t)  =  K(e2 ,  t) e 

is  given  by  equation  (25)  ,  the  connection 
a*  (e ,  t)  =  a  .  (e,  t-1e ) 

/  f  ^ 


(33) 


(34) 


now  follows  from  (31)  and  (33) .  It  follows  from  (34)  that 
if  3(e,  t)  takes  its  maximum  value  at  a  strain  of  1%  then 
so  does  o*(e,  t)  ,  i.e.  the  non-monotonicity  of  the  constant 
strain  rate  response  is  modelled  through  the  dependence  of  the 
kernel  on  strain,  the  value  of  1%  representing  a  failure 
strain  in  both  the  constant  strain-fate  and  constant  strain 
tests.  This  is  in  direct  contrast  to  the  linear  fluid  model 
(14) ,  where  the  mathematical  structure  forces  the  modelling  of 
the  non-monotonicity  of  the  constant  strain-rate  response  to 
be  done  through  the  time  dependence  of  the  kernel. 

Examination  of  equation  (31)  shows  that  the  kernel  K  can 
only  be  determined  by  the  constant  strain-rate  test  in  a 


region  of  the  (e ,  t) -plane  defined  by 


r  t  > 
max 


>  r  .  t 
mm 


(35) 


and 


s  e  , 
max  ' 


(36) 


where  e  is  the  maximum  attainable  experimental  strain 

max  r 

and  r  ..  and  r  .  represent  the  limits  within  which  the 
max  min  c 

strain  rate  falls. 

It  is  not  possible  to  evaluate  the  single  integral  relation 
of  Spring  and  Morland  (1983,  eq.  20,  p.189)  in  the  case  of 
the  constant  strain  test  because  this  model  will  not  admit 
a  strain  jump,  but  it  is  possible  to  examine  the  response  for 
the  strain  history 


r  (t) 


0  ,  t  <  0  , 

rt  ,  0  s  t  s  t^  , 

rtl '  tl  <  t  ' 


(37) 


where  the  time  t-^  is  small  and  r  is  constant.  The 
substitution  of  (37)  into  Spring  and  Morland’ s  equation  gives 

o(t)  =  T (r ,  t)  -  T ( r ,  t-t^  (38) 

for  t  >  t^  ,  where  T(r,  t)  is  the  stress  at  time  t  in 
a  uniaxial  compression  test  at  constant  strain-rate  r  . 


A  sketch  of  a  typical  stress  relaxation  curve  given  by 
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equation  (38)  is  presented  in  figure  2.  It  can  be  seen  from 
the  graph  of  o(t)  that  the  response  is  not  only  non-monotonic 
but  the  stress  actually  changes  sign  before  decaying  to  a 
constant  value  of  zero.  This  response  is  similar  to  the  constant 
strain  response  exhibited  by  the  "neo-Hookean"  fluid  model 

(14)  when  the  modelling  of  non-monotonic  constant  strain  rate 
behaviour  is  attempted.  Spring  and  Morland's  model  does, 
however,  give  a  good  account  of  the  variation  of  the  maximum 
strain-rate  but  due  to  the  kernel's  dependence  on  strain  rate 
rather  than  strain  precludes  the  modelling  of  non-monotonic 
constant  strain -rate  response  other  than  through  the  time 
dependence  of  the  kernel.  It  is  for  this  reason  that  the 
model  of  Spring  and  Morland  is  rejected  in  favour  of  the  present 
approach . 


5 .  Constant  Stress  Test.  I.  Approximate  Solution 

The  problem  of  the  constant  stress  stess  may  be  summarised 
in  the  following  question:  given  the  uniaxial  stress  history 


o  ( t) 


0  ,  t  <  0 

0  ,  t  £  0 


(39) 


• 

V.v 
•  »"■  - 
-y-\ 
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where  o  is  the  constant  stress,  what  form  does  the  strain 
e(t)  (t  >  0)  take?  An  approximate  solution  is  obtained  by 
first  addressing  the  more  general  question  of  finding  the 
strain  e(t)  (t  >  0)  generated  by  a  stress  field  which 

vanishes  for  t  <  0  but  is  arbitrary  for  t  >  0;  that  is, 

o (t)  =  0  ,  t  <  0  (40) 


...  ...  .  . .  _  »-  w... ..  —  .«L. 


Making  the  substitution 


u  =  -  efc (x) 


in  equation  (28),  noting  that  (40)  implies  that  e(t)  =  0 
for  t  <  0  ,  gives  for  t  >  0  , 

fe 

a(t)  =  (K[u2,t(e)-t  (e-u)  ]  +  2u2K  ,[u2  ,  t  (e ) -t  (e-u)  ]  }du  .  (42) 

>  1 


It  now  follows  from  (30)  that 


K(ez,t)  +  2e2K1(e2,t)  =  d^e, t  J'e) 


and  so  (42)  may  be  expressed  in  the  simpler  form 


t. 

■ 

o(t)  =  j  ojL:u,  t(c)  .ut(e.u|]  du  . 


Now  make  the  assumption  that  events  in  the  recent  past  dominate 


the  response  of  the  material  to  a  greater  extent  than  events 


in  the  distant  past.  Then 


t  (e)  -  t(e-u) 


e[  t(e) ] 


should  represent  a  reasonable  approximation  in  the  integral 
(44),  and  leads  to  the  approximate  relation 


Equation  (46)  is  a  first  order  ordinary  differential  equation 
for  the  strain  e  (t) .  It  may  also,  however,  be  regarded  as 
an  approximate  "equation  of  state",  replacing  the  integral 
relation  (23)  when  suitable  circumstances  apply. 

It  is  interesting  to  note  that,  when  the  stress  history 
(39)  is  employed  in  (46)  ,  the  constant  stress  data  is  given 
by  the  one-parameter  family  of  curves  in  the  (e ,  e)  plane 
represented  by 


o  (e  ,  e  ) 


=  a 


(47) 


a  being  the  parameter.  Furthermore,  since  the  constant 
strain  rate  data  is  represented  by  the  family 


a  =  a  (e ,  r) 


(48) 


of  curves  in  the  (a,  e)  plane,  it  can  be  verified  that  the 
conjectures  of  Mellor  and  Cole  (1982,  Table  1,  p.206),  connecting 
the  constant  strain  rate  and  constant  stress  data,  are  borne 
out  by  (47)  and  (48).  Furthermore,  the  relationship  between 
the  constant  stress  and  constant  strain  rate  data  embodied 
in  (47)  and  (48)  has  recently  been  verified  experimentally  by 
Mellor  and  Cole  (1983)  to  a  reasonable  degree  of  accuracy. 

We  are  thu~  led  to  the  conclusion  that  the  constant  strain 


rate  and  constant  stress  data  represent  essentially  the  same 
information. 


Equation  (47)  is,  however,  only  an  approximate  solution  to 


the  problem  of  the  constant  stress  test  and  a  better 
correlation  between  the  constant  stress  data  and  the  constant 
strain-rate  data  should  be  obtained  from  a  numerical  solution 
of  the  integral  equation  associated  with  the  constant 
stress  test. 


6.  The  Structure  of  the  Kernel 


In  this  section  the  mathematical  structure  of  the  kernel 
K(e2,t)  is  examined  with  the  aid  of  the  assumption  that 
the  constant  strain-rate  data  can  be  expressed  in  a  separable 


form 


5{e,  r)  =  s(e)amax(T>  ' 


This  assumption  is  consistent  with  the  experimental 

observation  that  the  stress-strain  curves  for  different 

constant  strain-rates  are  non-intersecting.  The  function 

S(e)  is  referred  to  as  the  shape  function  and  is  sketched 

in  figure  3.  a  (r)  is  the  maximum  stress,  which  is 

max 

sketched  in  figure  4.  It  should  be  noted  that  no  attempt 
has  been  made  to  include  the  experimentally  observed  fine 
structure  to  the  left  of  the  peak  in  the  shape  function  because 
this  is  believed  to  be  due  to  internal  cracking  (Mellor  and 
Cole,  1982),  which  is  outside  the  scope  of  this  present  model. 

When  (49)  is  substituted  into  (34)  and  the  result  is 
integrated  by  parts  the  relation 


a*(e,t)  =  S  (e)  a _ (t  e) 


-  [£  S(u)o _ '  (t"1u)t"1du 


(50) 


is  obtained.  When  the  substitution  tv  =  u  is  made  in  the 


integral  on  the  right-hand  side  of  (50)  it  is  apparent  that 

o* (e,0)  =  omS(e)  .  (51) 

By  definition,  o*(e,  0)  is  the  instantaneous  elastic 

response  and,  apart  from  the  constant  multiplier 

a  =  lim  o  (r)  ,  it  is  the  same  as  the  shape  function  so 
00  r-*-oo  max 

that  reference  may  again  be  made  to  figure  3.  It  is  now 
possible  to  sketch  the  initial  value  K(ez,0)  of  the  kernel 
and  this  is  given  in  figure  5. 

It  may  also  be  inferred  from  the  result 

a*,1(e,t)  =  S' (e)omax(t'1e)  ,  (52) 

which  is  a  consequence  of  (50),  that  a*,^(e,t)  (t  >  0)  is 
of  a  similar  nature,  viewed  as  a  function  of  e  ,  to 
o*,^(e,0),  with  the  exception  that  the  former  vanishes  at 
z  -  0  .  Therefore  it  follows  that  K(e2,t),  which  is  given 
by  (33),  has  a  similar  shape  to  the  ini-tial  kernel  K(e2,0) 
shown  in  figure  5 . 
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